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Abstract 

Entropy production from scalar decay in the era of low temperatures after neutrino decoupling 
will change the ratio of the relic neutrino temperature to the CMB temperature, and, hence, the 
value of Nf,ff, the effective number of neutrino species. Such scalar decay is relevant to reheating 
after thermal inflation, proposed to dilute massive particles, like the moduli and the gravitino, 
featuring in supersymmetric and string theories. The effect of such entropy production on the relic 
neutrino temperature ratio is calculated in a semi-analytic manner, and a recent lower bound on 
this ratio, obtained from the WMAP satellite and 2dF galaxy data, is used to set a lower bound 
of ~ 1.5 X 10"^'^ Gev on the scalar decay constant, corresponding to a reheating temperature of 
about 3.3 Mev. 
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I. INTRODUCTION 



The question of neutrino equilibration and subsequent decoupling in the early universe 
has become increasingly important over the past few years with the advent of precision 
cosmology 

The canonical textbook result for neutrino decoupling is that neutrinos decouple prior to 
electron-positron annihilation, leading to a final neutrino temperature which is related to 
the photon temperature by r = T^q/Tq = (4/11)^/'^, where T^q^Tq, are, respectively, the relic 
neutrino and CMBR temperatures. 

There are quite a few factors which may cause slight departure from the standard value of 
r = (4/11)^/'^. First, the neutrinos, too, are slightly heated during the e~^e~ annihilations jj], 
leading to an overall increase in neutrino energy density of roughly 1%. Finite temperature 
QED effects lead to an addMonal =UgM heating of the neutrinos k Howevet, all ot these 
results are within the standard model and lead to an overall increase in neutrino energy 
density of a little more than 1%, an effect which should be included in practical calculations. 

At present this small effect is not detectable, but it has been estimated that future high- 

aecision measurements of the CMB anisotropy could reach the required level of sensitivity 

In the present paper we discuss an additional factor which has effect on the value of 
r. Various long-lived, massive fields like the gravitino, the Polonyi, the moduli, and the 
dilaton. wh,ch figure m supersyutmetdc and stnug theory models fl. pose cosmological 

problems, because their decay must affect rj and nucleosynthesis [5|]. To dilute them away, 
the proposal of thermal inflation has been mooted. A scalar fleld, the flaton, is used to 
generate inflation at late times, typically at temperatures of about 10^ Gev. The inflation 
stops when the temperature falls to the flaton mass ~ 10^ Gev. Such a particle will go on 
decaying into the era of Mev-scale temperatures, and affect nucleosynthesis and the CMBR. 

The effect on neutrino decoupling and nucleosynthesis has been studied 0, 0, ll^ • The 
studies do not depend materially on the details of flaton phenomenology, and apply to the 
Mev-scale decay of any scalar 0. 

In this manuscript, we are trying to see how the outpouring of entropy, as decays, heats 
up the e~,e"'",7 plasma, and affects the ratio r. The assumption of current models is that 
the (p should not decay into neutrinos directly [7] . So, the decay of (p operates in a direction 
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opposite to neutrino heating, and may offset conclusions drawn on the basis of neutrino 
heating. The effect of changing the lower bound on the scalar decay constant F should be 
an interesting input in current calculations of nuclear abundances and CMBR anisotropies. 

Ref. j?! deals directly with the neutrino distribution function and proceeds by solving 
the Boltzmann equation numerically. We try, on the other hand, to adhere as far as possible 
to the macroscopic entities like temperature and entropy so as to proceed analytically and 
keep the physical processes transparent. 

Of course, this means that our calculation does not reach the level of precision of a 
numerical solution to the Boltzmann equation. However, it does provide a reasonable bound 
on the involved parameters. 

The plan of the paper is as follows. Section 1 is this Introduction. In section 2, the 
entropy production due to scalar decay, since neutrino decoupling, is estimated, and the 
basic relation whereby it affects the relic neutrino temperature ratio is set out. In section 3, 
the values of this ratio for different values of the scalar decay constant are calculated. Section 
4 uses a lower bound on the relic neutrino temperature ratio from the WMAP satellite data, 
the 2dF galaxy survey and other related data to find a lower bound on the scalar decay 
constant, and discusses the results. 



II. CALCULATION OF THE ENTROPY PRODUCTION 

Kawasaki et al [3] find that the electron and muon/taon neutrino distribution functions 
show a deficit from the thermalised F.D. distribution as the reheating temperature for the 
scalar (and, hence, its decay constant) falls. This leads to a decrease in the Ve and v^^r 
energy densities with concomitant effect on the weak interaction rates and freeze-out times. 
The threshold reheating temperature Tr is around 7 Mev, below which the authors find the 
effective number of neutrino types N^ff falls below the value 3. 

So, even before decoupling, the neutrinos cannot be assigned the photon temperature. 
Can they be assigned a temperature at all, in particular a lower temperature which will 
approximately reproduce the decreased distribution? Strictly speaking, the shape of the 
distribution does not permit this. Ref. j3] defines = [27T'^n^/ {3({3))]''^^^^ and Re = 
{pu / nu) / {3.151Ti^) , and finds that Re takes values 1.00, 1.03, 1.50 for Tr values 10,3, 1 Mev, 
respectively, while a F.D. distribution should give Re = 1- On the other hand, these 



3 



results indicate that unless the decay constant F is much smaller than that required to save 
standard BBN, the assumption of a neutrino temperature T^, a little smaller than the photon 
temperature T, even before decoupling, is not a bad one, especially if one is interested only in 
locating the parameter region where F does not affect the standard picture. After decoupling, 
the neutrinos are no longer affected by the effects of scalar decay, and the momentum of 
the neutrinos, assumed here to be massless, redshifts as 1/a, so that one can consider the 
usual relic neutrino temperature-like parameter which redshifts as 1/a from the neutrino 
decoupling temperature Tj^j. 

However, the main difference between T^, and T must arise from the e~, annihilations, 
after the decoupling era, and we are interested in seeing how this is affected by scalar decay. 
In what follows, we will neglect the difference between T^i and Tj, just at the epoch of 
decoupling t — U, while recognising that the decoupling temperature remains uncertain to 
some extent. To take this into account, we will consider a decoupling temperature range of 
1 - 3 Mev. 

We have to deal with four epochs, i is an epoch just after electron neutrino, v, decoupling, 
when we take T^i = Ti, and 

g*s, « 9* = 11/2. (1) 

/ is an epoch after e~ , annihilation, which we can take to be the present, with neutrino 
temperature Ti,q and photon temperature To, and 

9*Sf = = 2. (2) 

a is the epoch of e~,e"'" annihilation, when the photon temperature is T^, and, although in 
this epoch, the effective number of degrees of freedom is actually changing, we will take 

g*sa^9:^9sr^9:^lll'^- (3) 

A basic assumption is that at some earlier era, dominated the energy density. We 
will need a fiducial era for 0, when this 0-domination of the universe ends and radiation 
domination begins. This, we assume, happens sometime before decoupling of the three 
families of neutrinos, at a temperature Tg, with g*g^ = 43/4. 

The basic relation is 

gU'^^'/^h)a\Tf + = ^J/27rV45)ajT|. (4) 
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AS is the entropy poured into the plasma between Tj and Tf, in accordance with the 
assumption that 0-decay into neutrinos is not allowed. It is to be calculated from 

dS = -d{a^p^)/T, (5) 

P0 being the scalar energy density at radiation temperature T, radiation including all mass- 
less particles in equilibrium with photons. If pr is the radiation density, temperature is 
defined from 

PR = {ny30)g*T\ (6) 
the equation for the evolution of the 



Taking the scalar decay rate constant to be F ^| 
scalar energy density is 



d 

—p^ + 3Hp^ = -Tp^. (7) 
Defining $ = a^p,^ and R = a^pn, where a is the scale factor, becomes 

$ = -r$, (8) 

with the solution [3] 

$ = ^^e-'^^*-*^), (9) 

where is the value of $ at t = tg. So, 

S = {r/T)<l>Ee-^^'-'^\ (10) 

In (jlOj) . the exponential will dominate at the epochs i and /, and in between, when pn 3> 
Pif,. So, in the pre-exponential, we may approximate T from © using the full radiation 
domination equations H = l/(2t) and = 8ttpr/{3Mpi). This gives the usual relation 

1/T = aVi, with = 2.7215 x 10'^^g*^2Gev~\ (11) 

Introducing the variable 

y = Tt = T/ia'T'), (12) 

we get 

rf 

AS = r$se^*^ / aVie-^'dt 

J i 

= {l/Vf)^Ee^'^ J^^ a^e-ydy 

J i 

+ V^e-ydy)l (13) 



where a, is the value of a with g* = g* and af/ai = We have assumed g* to be g* 

from ti to ta and gj from to tj. 



A. Estimate of 

Because there are as yet no firm phenomenological values for ^e, we have to make an 
estimate. In the expression for S given in ()1U|) . the exponential dominates at the epochs we 
are interested in. So, we estimate in the pre-exponential just to an order of magnitude, 
defining the epoch E by taking = at t = Ie- Just for estimating in the pre- 
exponential, we use the crude approximation a oc 1/T, although, actually, in the presence of 
entropy generation there is departure from this typ e of evolution, and, in fact, in the regime 
of full $ domination, when p^ ^ pR, T oc a~i Iii the regime of pr ^ p^, of course, 

we expect a closer fit to a oc 1/T. 

With this approximation, in the pre-exponential, 
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= O'ePre 

= {n'/30)g*Ea%Tl. 

= {7r'/30)g*EafT^TE. (14) 



B. Estimate of Te 

Te is the temperature when ^-domination passes into /^-domination, i.e. when passes 
from -^^7-P(j, to -^^Pb- So, we try to find approximations for closer to the epoch p^ = p^, 
and on either side of it. We first consider the era before this, i.e. the era of incomplete $- 
domination. From 

d d 
^[a'(p0 + Pi?)]+PK^a' = O, 

and (IZI),®, one obtains 

R = ar$. (15) 
The Friedmann equation can be written as 



In the era of incomplete ^-domination, the term £^ on the RHS of p6|) is a correction term, 
and may be evaluated to the approximation 

dR 

R- Ri = -^\ai{a- ai), (17) 

where tj refers to some initial epoch such that a ^ a/. Also, if it is supposed that the 
scalar decay produces sufficiently copious radiation, Rj <ti R. As a correction term is being 
dealt with, these approximations should not cause much deviation from the actual evolution. 
Then, for t sufficiently later than tj, but within the regime under consideration, one may 
write, in the correction term on the RHS of (fTT)|) . 

r 

R ^ — $a, (18) 
H 

using ()15|1 and p7|l . As we are hoping to find only an order of magnitude estimate of a 
lower bound on F, (fTHj) is not a bad approximation. Thus, if we go even to very early times 
in this era, when $ = constant 13], integration of ()15|) leads to 

2 r 

using the approximations a ~ ti, H = ^. 

So, introducing a new evolution variable x = r/H, we write, for use in the correction term 
on the RHS of (HH), 

x = T/H^ R/^a = pn/p^. (19) 

Next, we consider the era of interest to us when $a <^ i?, such that ^a/R cannot be 
neglected, but its higher powers can. In this era of incomplete radiation domination, the 
Friedmann equation is put in the form 

If, well into this epoch, the correction term ^a/R on the RHS of (PUI) is neglected, the full 
radiation domination relations are found: 

H = — , and 

a = At^, (21) 

A being a constant. 
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(jS} has, as solution, a falling exponential in t, viz. $ ~ e~^*. Instead of taking the falling 
exponential in t directly, a suitable approximation to the correction term on the RHS of (j^Uj) 
is first worked out. Let to be a sufficiently late epoch, when $ = $o ~ 0. Then, for use only 
in the correction term on the RHS of (j2(Jj) . one takes 

*-4„ = 4(-)-4(_, = -^|..(-__,. 

Neglecting $o, I /to compared to l/t, respectively, an approximation 

$ ^ P (22) 

will be used only in the correction term on the RHS of ()2U|) . i.e. in the correction term, 
the falling exponential will be approximated by a rectangular hyperbola. B is a constant. 
A similar approximation is considered for R. It ought to be mentioned that R refers to the 
total radiation present, and not only to that produced by decay. However, the change in R 
is due to decay and consequent entropy production. In the absence of this decay, R = 0. 

Using ()21|) and ()22|) in (fT3j) . and, integrating, one obtains, for use only in the correction 
term on the RHS of 

R-Re ^2ABT{t^ -4), 

an approximation which corresponds to (fT7j) . because of a ~ At^. If t^; is sufficiently 
early compared to t, though within the regime under consideration, and there is sufficiently 
copious radiation production since Ie-, it is sufficient to take 

R ^ 2ABn^ 

in the correction term on the RHS of ((201) • This relation, together with (^1]) and (j22I), are 
now used to give, in the correction term on the RHS of (j2Up . 

r 



once again, as in (fT^ . 

Introducing the variable x in ()20j) . we get ()20p in the form 

8ttR ,^ 1. 



(23) 



a;2 SMV^^ + x^- ^^^^ 



This is a good equation for incomplete radiation domination when x ^ 1. However, we will 
approximate Te in i.e. in the pre-exponential of (jTU)), by putting x = xe ^ ^ for t = tE 
in (j211). 

Now, (fTUI) and signify that our approximations are equivalent to taking p$ = 
when r = if, to an order of magnitude, in the correction term in H^, and this is our way 
of bypassing lack of knowledge about the initial value of $. This approximation has been 
explained in the preceding parts of this section. It differs from a common approach to the 
problem where decay is supposed to occur at t = l/F and $ is put equal to the scalar mass 
at this epoch. Here, we allow the scalar to decay over time, but the price of bypassing 
knowledge about scalar energy density or mass is paid by the approximation inherent in 
(fTI?|) and (f2H|) and our way of estimating Te- The result for Te is 

TE^2-^y/f/aE. (25) 

If we use fll9|) in the incomplete 0-domination case when x <^ 1, equation (fTBj) becomes 

-x). (26) 



x2 3M|,;a3 



This gives the same value for Te if this equation is extrapolated to x = xe ^ ^ ioi t = Ie- 
Putting the value of Te from (^3)) in (HU, (fTSj) becomes 



4.1,/-/ 

+ ( 

Now, in dH), (EZj) 



AS = 2^{'K^/'iQ)g],a^,Tfei{a,/aE)[j^ ^e-^dy 

+ V^e-ydy)l (27) 



a:^f = a}T^,, (28) 

because the temperature of the decoupled neutrinos red-shifts as 1/a. Using and (|77j). 
we get 

+ 2A31r'[J^^ ^e-ydy + ^e-ydy)] = A/11, (29) 

where we have put r = T^q/Tq, and taken (XijaE = (di/gE)^ = (22/43)3. 



III. RESULTS 



From (|I2I),(IIII),(II1),(|1 and ©, we find 



_ 1.5667510 _ 2.60ro _ 1.5667510 , . 

~ {Ti/Mevy^^ ~ {Tf/Mevy'^" ~ {Ta/Mevf ^ ' 
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Ti = IMev 


Ti = 2Mev 


Ti = 3Mev 


To 


r 


To 


r 


To 


r 


1.0 


0.587 


4 


0.587 


10 


0.599 


1.2 


0.608 


5 


0.612 


12 


0.620 


1.4 


0.626 


6 


0.635 


13 


0.630 


1.5 


0.635 


7 


0.654 


14 


0.639 


1.6 


0.643 


8 


0.669 


15 


0.648 


1.8 


0.657 


10 


0.690 


16 


0.656 


2.0 


0.669 


12 


0.702 


20 


0.680 


3.0 


0.702 


14 


0.708 


24 


0.695 


4.0 


0.711 


20 


0.713 


28 


0.703 


10 


0.714 


25 


0.714 


32 


0.708 


25 


0.714 


50 


0.714 


50 


0.713 










75 


0.714 










100 


0.714 



TABLE I: Relic Neutrino Temperature Ratio r for different values of the Scalar Decay Parameter 
To = r/{10-^^Gev) 

Ti is a few Mev. Here, we calculate results for Tj = 1,2,3 Mev. T/ being << leV, yf can 
be put ~ oo because of the nature of the incomplete Gamma function r(1.5, x), where 

ry2 

/ ^e-ydy = T{1.5,y,)-T{1.5,y2). (31) 
Jyi 

Also, the properties of the incomplete Gamma function r(1.5,x) indicate that the results 
should be insensitive to the precise value of in the range 0.3 < Ta/Mev < 0.5. We have 
checked this numerically for Ta/Mev = 0.3,0.4,0.5. The values of the ratio r for different 
values of Tq are displayed in Table 1. 

IV. CONCLUSIONS 

A semi-analytical method for calculating the change in the relic neutrino temperature 
ratio and, hence, in the effective number of neutrino species, resulting from reheating at 
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very low temperatures, has been presented. 

While this calculation obviously does not have the accuracy of a full numerical solution of 
the Boltzmann equation, it does have the merit of being transparent and easily reproducible. 

In Ref. , a lower bound on the reheating temperature was calculated from a considera- 
tion of the impact of incomplete neutrino equilibration on big bang nucleosynthesis. While 
this bound is powerful, it has the problem that it is flavour dependent in the sense that 
incomplete electron neutrino equilibration has a direct impact on the n — p conversion rate. 

Here we use, instead, a bound which relies on energy density only. Recent results [l^ 
show that an overall best fit for the WMAP T and TE data, combined with the Wang et 
al compilation of CMB data [l^, the 2dFGRS data, the HST key project data on Hq, and 
the SNI-a data on fi^, give bounds on A^e// : 

1.9 < Neff < 7.0 (95% confidence ). (32) 

This corresponds to a lower bound on r : 

r > 0.637. (33) 

Table I shows that the change in r due to F is sensitive to Tj, i.e. the neutrino decou- 
pling temperature. The lower bound on F for a lower bound of 0.637 on the relic neutrino 
temperature ratio, corresponding to (jH^ . is found to be between 1 and 2 x 10~^^Gef , 5 and 
7 X lO-^^'Get;, 12 and 15 x lO-^^'Gew, respectively, for = 1, 2, 3 Mev. We conclude that to 
ensure that the relic neutrino temperature ratio r = T^q/Tq does not exceed the lower bound 
of 0.637 1^, the scalar decay constant F must be greater than about 15 x 10~^^G'ef . Taking 
the definition of the reheating temperature in Q, this corresponds to a reheating tempera- 
ture of 3.3 MeV, which is comparable to the result found in Ref. 0. In their calculation, 
N(.ff = 1.9 corresponds roughly to Tr ~ 2.2 MeV. 

This means that even though the approximations we use are rough, the end result is very 
similar to what is found from the full numerical solution. The reason is that we are only 
interested in energy density, not in the underlying neutrino distribution function. On the 
other hand, for BBN purposes the distribution function of Vf. is very important, because 
high energy neutrinos have more weight in the n ~ p conversion processes. It should also 
be noted that the lower bound on Tr of 0.7 MeV found in Ref. Q is not based on N^ff 
alone, but rather on a detailed study of primordial abundances. This clearly illustrates 
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the fact that BBN bounds are highly flavour sensitive, and shows that the much simpler 
energy density bound from CMB and large scale structure leads to a stronger bound on the 
reheating temperature. 
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